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ABSTRACT 

The extent of nucleosynthesis in models of type I X-ray bursts and the associated impact on the 
energy released in these explosive events are sensitive to nuclear masses and reaction rates around the 
64 Ge waiting point. Using the well known mass of 64 Ge, the recently measured 65 As mass, and large- 
scale shell model calculations, we have determined new thermonuclear rates of the 64 Ge(p, 7 ) 65 As and 
6o As(p, 7 ) 66 Se reactions with reliable uncertainties. The new reaction rates differ significantly from 
previously published rates. Using the new data we analyze the impact of the new rates and the 
remaining nuclear physics uncertainties on the 64 Ge waiting point in a number of representative one- 
zone X-ray burst models. We find that in contrast to previous work, when all relevant uncertainties 
are considered, a strong 64 Ge rp-process waiting point cannot be ruled out. The nuclear physics 
uncertainties strongly affect X-ray burst model predictions of the synthesis of 64 Zn, the synthesis 
of nuclei beyond A = 64, energy generation, and burst light curve. We also identify key nuclear 
uncertainties that need to be addressed to determine the role of the 64 Ge waiting point in X-ray bursts. 

These include the remaining uncertainty in the 65 As mass, the uncertainty of the 66 Se mass, and the 
remaining uncertainty in the 65 As(p, 7 ) 66 Se reaction rate, which mainly originates from uncertain 
resonance energies. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: neutron — X-rays: bursts 


1. INTRODUCTION 

A type I X-ray burst (XRB) arises from a thermonu¬ 
clear runaway in the accreted envelope of a neutron star 
in a close binary star system (for reviews, see, e.g., Lewin 
et al. (1993); Schatz et al. (1998); Strohmayer & Bild- 
sten (2006); Parikh et al. (2013)). Roughly 100 bursting 
systems have been discovered to date, with light curves 
exhibiting peak luminosities of Apeak ~ 10 4 -10 5 Lq and 
timescales of 10-100 s. During an XRB, models predict 
that a H/He-rich accreted envelope may become strongly 
enriched in heavier nuclei though the op-process and 
the rp -process (Wallace & Woosley 1981; Schatz et al. 
1998). These two processes involve o-particle-induced or 
proton-capture reactions on stable and radioactive nu¬ 
clei, interrupted by occasional /3-decays. When the rp- 
process approaches the proton dripline, successive cap¬ 
ture of protons by nuclei is inhibited by a strong reverse 
photodisintegration reaction rate. The competition be¬ 
tween the rate of proton capture and the rate of /3-decay 
at these “waiting points” (e.g., 60 Zn, 64 Ge, and 68 Se) de¬ 
termines the extent of the synthesis of heavier mass nu¬ 
clei during the burst (Schatz et al. 1998). Peak temper¬ 
atures during the thermonuclear runaway may approach 
or exceed 1 GK, resulting in the synthesis of nuclei up to 
mass A « 100 (Schatz et al. 2001; Elomaa et al. 2009). 
Model predictions depend, however, on astrophysical pa- 
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rameters such as accretion rate, the composition of the 
accreted material, and the neutron star surface gravity, 
as well as on nuclear physics quantities such as nuclear 
masses and reaction rates. 

The 64 Ge(p, 7 ) 65 As and 65 As(p, 7 ) 66 Se reactions have 
been demonstrated to have a significant impact on nu¬ 
cleosynthesis during XRBs. (See Parikh et al. (2014) for 
a recent review of the impact of nuclear physics uncer¬ 
tainties on predicted yields and light curves from XRB 
models.) Direct measurements of these reactions at the 
relevant energies in XRBs are not yet possible due to 
the lack of sufficiently intense radioactive 64 Ge and 65 As 
beams. Moreover, due to the unknown mass of 66 Se and 
the lack of nuclear structure information for states within 
r; 1-2 MeV of the 64 Ge+p and the (theoretical) 65 As±p 
energy thresholds in 65 As and 66 Se, respectively, it is not 
possible to estimate rates for these reactions based on 
experimental nuclear structure data. As a result, XRB 
models use 64 Ge(p, 7 ) 65 As and 65 As(p, 7 ) 66 Se thermonu¬ 
clear rates derived from theoretical calculations. Using 
such models it has been demonstrated that varying the 
65 As(p, 7 ) 66 Se rate by a factor of ten at the relevant tem¬ 
peratures affects the calculated abundances of nuclei be¬ 
tween A « 65-100 by factors as large as about 5 (Parikh 
et al. 2008). For 64 Ge(p, 7 ) 65 As, models have illustrated 
the importance of the Q-value (or proton separation en¬ 
ergy S p ) adopted for this reaction, with variations by 
±300 keV affecting final calculated abundances between 
A rs 65-100 by factors as large as about 5 (Parikh et 
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al. 2008, 2009). In addition, the effective rp- process 
lifetime of the waiting-point nucleus 64 Ge was investi¬ 
gated by Schatz (2006) based on the estimated proton 
separation energies of S' p ( 65 As)=-0.36±0.15 MeV and 
S’ p ( 66 Se)=2.43±0.18 MeV, derived from Coulomb mass 
shift calculations (Brown et al. 2002). It was found that 
the effective lifetime of 64 Ge for a given temperature and 
proton density is mainly determined by the S p values of 
65 As and 66 Se and the proton capture rate on 65 As. 

Recently, precise mass measurements of nuclei along 
the rp-process path have become available. The mass of 
64 Ge has been measured at the Canadian Penning Trap 
at Argonne National Laboratory (Clark et al. 2007) and 
the LEBIT Penning Trap facility at Michigan State Uni¬ 
versity (Schury et al. 2007). More recently, the mass 
of 65 As has been measured at the HIRFL-CSR (Cooler- 
Storage Ring at the Heavy Ion Research Facility in 
Lanzhou) (Xia et al. 2002) using IMS (Isochronous Mass 
Spectrometry). The measurements can be combined to 
obtain an experimental proton separation energy for 65 As 
of S p = —90 ± 85 keV (Tu et al. 2011), where the un¬ 
certainty is dominated by the uncertainty in the 65 As 
mass. The mass of 66 Se is not known experimentally. 
The extrapolated value predicted by AME2012 results 
in S'p( 66 Se) = 1720 ± 310 keV. With the new mass of 
65 As, X-ray burst model calculations (Tu et al. 2011) 
suggested that 64 Ge may not be a significant rp-process 
waiting point, contrary to previous expectations (Schatz 
et al. 1998; Woosley et al. 2004; Fisker et al. 2008; Parikh 
et al. 2009; Jose et al. 2010). We revisit this question here 
using our new nuclear reaction rates. 

Thermonuclear 64 Ge(p,q) and 65 As(p,q) reaction rates 
were first estimated by Van Wormer et al. (1994) based 
entirely on the properties of the mirror nuclei 65 Ge and 
66 Ge, respectively. S p values of 65 As and 66 Se were es¬ 
timated to be 0.169 MeV and 1.909 MeV, respectively. 
Later on, both rates have been calculated (Rauscher & 
Thielemann 2000) with the statistical Hauser-Feshbach 
formalism (NON-SMOKER (Rauscher & Thielemann 
1998)) using the masses of 65 As and 66 Se predicted by 
the finite-range droplet (FRDM) (Moller et al. 1995) and 
ETSFIQ (Pearson et al. 1996) mass models. Recently, 
the statistical model calculations have been updated us¬ 
ing new predictions for the 65 As and 66 Se proton separa¬ 
tion energies (see JINA REACLIB 1 (Cyburt et al. 2010)). 
The predicted rates differ from one another by up to sev¬ 
eral orders of magnitude over typical XRB temperatures. 
Moreover, the reliability of statistical model calculations 
for these rates is questionable due to the low compound 
nucleus level densities, especially for 64 Ge(p,q), but also 
for 65 As(p,q). 

In this work we refer to previously available rates us¬ 
ing the nomenclature adopted in the JINA REACLIB 
database. The laur rate refers to the rate estimated 
by Van Wormer et al. (1994); the rath rate was calculated 
by Rauscher & Thielemann (2000). The rath, thra, rpsm 
rates are the statistical-model calculations with FRDM, 
ETSFIQ, as well as Audi & Wapstra (1995) estimated 
masses, respectively. The recent ths8 rate is from Cy¬ 
burt et al. ( 2010 ). 


Here we determine new thermonuclear 64 Ge(p,q) 65 As 
and 65 As(p,q) 66 Se reaction rates using the updated S p 
values of 65 As and 66 Se together with new nuclear struc¬ 
ture information from large-scale shell-model calcula¬ 
tions. Using the new data we fully characterize the 
nuclear physics uncertainties that affect the rp-process 
through 64 Ge and reexamine the question of the 64 Ge 
waiting point. 

2. REACTION RATE CALCULATIONS 

The total thermonuclear proton capture reaction rate 
consists of the sum of resonant- and direct-capture (DC) 
on ground state and thermally excited states in the tar¬ 
get nucleus, weighted with their individual population 
factors (Fowler et al. 1964; Rolfs & Rodney 1988). It 
can be calculated by the following equation: 

(n j. _i_ i \ p -Ei/kT 

n a(<jv) = + N A (<rv)j3c) ^2J n + ly-En/kT 

( 1 ) 

with the parameters defined by Schatz et al. (2005). 


2 . 1 . Resonant rates 

For isolated narrow resonances, the resonant reac¬ 
tion rate for capture on a nucleus in an initial state 
i, Na^v)^, can be calculated as a sum over all rele¬ 
vant compound nucleus states j above the proton thresh¬ 
old (Rolfs & Rodney 1988; Iliadis 2007). It can be ex¬ 
pressed by the following equation (Schatz et al. 2005): 


N A (<rvy r = 1.54 x 10 n (pT 9 )~ 3/2 

v- ( 11.605FU 

x 0 ^ 7 ^-exp I - 

j ^ 


To 


^ [cm 3 s 4 mol 4 ], 
( 2 ) 


where the resonance energy in the center-of-mass system, 
Eij = Ej — S p — Ei, is calculated from the excitation 
energies of the initial Ei and compound nucleus Ej state. 
For the ground-state capture, the resonance energy is 
represented by U* = E° x — S p . T 9 is the temperature 
in Giga Kelvin (GK) and p is the reduced mass of the 
entrance channel in atomic mass units (p = At/{ 1 + 
At), with At the target mass number). In Eq. 2, the 
resonance energy and strength are in units of MeV. The 
resonance strength wq is defined by 


2 Jj + 1 T*J x VI 
" 7 «-2(2J, + l) q ot>l 


( 3 ) 


where A is the target spin and Jj, T p 7 , T-L and T:( otal 
are spin, proton decay width, q-decay width, and total 
width of the compound nucleus state j, respectively. The 
total width is given by T 4 otal = W + because 

other decay channels are closed (Audi et al. 2012) in the 
excitation energy range considered in this work. 

The proton width can be estimated by the following 
equation, 


T p = ^C 2 S(nQ) T sp (nlj), 

nlj 


1 http://groups.nscl.msu.edu/jina/reaclib/db 


( 4 ) 
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TABLE 1 

Properties of 65 As for the ground-state capture utilized in the present 64 Ge(p,7) 65 As resonant rate calculation. The 
ENERGY LEVELS IN THE MIRROR NUCLEUS 65 Ge ARE LISTED IN THE 4 th COLUMN FOR COMPARISON. 



ET V [ a] 

E x [MeV] 
£‘ heo [b] 

S,( 65 Ge)[c] 

E r (MeV)[d] 

nlj 

C 2 S 

T P [eV] 

r 7 [eV] 

017 [eV] 

3/27 

0.000 

0.000 

0.000 

0.090 

2^3/2 

0.196 

1.19 x 10 " 34 

0.00 

0.00 

5/27 [e] 

0.187(3) 

0.103 

0.111 

0.277 

I/ 5/2 

0.533 

8.19 x 10 - 17 

5.11x 10 -7 

2.46xl0 -16 

5/27 


0.501 

0.605 

0.591 

I/ 5/2 

0.010 

3.76 x lO " 10 

4.73x 10 -5 

1.13x 10 -9 

5/27 


0.863 


0.953 

I/ 5/2 

0.014 

1.64 x 10 “ 6 

4.24x 10 —4 

4.89x 10 -6 

7/27 


0.947 

0.890 

1.037 

I/ 7/2 

0.013 

1.28 x 10“ 5 

2.43x 10 -4 

4.88X10 " 5 

7/27 [f] 


1.070 

1.155 

1.160 

I/ 7/2 

0.002 

8.50 x lO " 6 

2 . 11 x 10 -4 

3.27X10" 5 


a] measured by Obertelli et al. (2011); 

b] calculated by the present shell model; 

c] compiled by Browne & Tuli (2010); 

d] calculated by E r = E x — S p with S p = —0.09 MeV (Tu et al. 2011); 

e] calculated E r and F p based on the experimental value of E x = 0.187 MeV for this state; 

f] negligible contribution to the rate for temperatures up to 2 GK. 


where C 2 S(nlj ) denotes a proton-transfer spectroscopic 
factor, while r sp is a single-proton width for capture 
of a proton on an ( nlj ) quantum orbital. The r sp are 
obtained from proton scattering cross sections calcu¬ 
lated with a Woods-Saxon potential (Richter et al. 2011; 
Brown 2014). Alternatively, the proton partial widths 
may also be calculated by the following expression (Van 
Wormer et al. 1994; Herndl et al. 1995), 


The 7 widths, T t , have been calculated from the 
electromagnetic reduced transition probabilities B(Ji —> 
Jf,L ), which carry the nuclear structure information of 
the resonant states and the final bound states (Brussaard 
& Glaudemans 1977). The reduced transition rates are 
computed within the shell model. Most of the transi¬ 
tions in this work are of Ml and E 2 types. The relations 
are (Herndl et al. 1995): 


r p = jjpPt(E)C 2 S. (5) 

Here, R = ro x (1 + At) 1 ^ fm (with ro = 1.25 fm) is the 
nuclear channel radius. The Coulomb penetration factor 
Pt is given by 


Pt{E) 


kR 

F 2 {E) + G 2 {EY 


( 6 ) 


where k = y/2jlE/h is the wave number with energy E in 
the center-of-mass (c.m.) system; Ft and Gt are the reg¬ 
ular and irregular Coulomb functions, respectively. The 
proton widths given by these two methods (i.e., by Eqs. 4 
and 5) agree well with each other, with a maximum dif¬ 
ference of about 35%. 

The key ingredients necessary to estimate the reso¬ 
nant 64 Ge(p, 7 ) and 65 As(p, 7 ) rates are energy levels 
in 65 As and 66 Se, proton transfer spectroscopic fac¬ 
tors, and proton and gamma-ray partial widths. For 
65 As, only a single level has been observed at E x = 
187(3) keV (Obertelli et al. 2011). For 66 Se, one level 
has been confirmed at E x = 929(7) keV, and indications 
for two other levels at 2064(3) keV and 3520(4) keV have 
been reported, with tentative assignments of (4 + ) and 
( 6 + ), respectively (Obertelli et al. 2011; Ruotsalainen 
et al. 2013). There are no more experimental data 
available for these two nuclei. In this work, we have 
calculated the energy levels, spectroscopic factors and 
gamma widths within the framework of a large-scale shell 
model, without truncation, using the shell-model code 
NuShellX@MSU (Brown & Rae 2014). The effective in¬ 
teraction GXPFla (Honma et al. 2004, 2005) has been 
utilized for these two p/-shell nuclei. 


r E2 [eV] = 8.13 x 10” 7 ^[MeV]H(F2)[e 2 fm 4 ], (7) 

and 

r M i[eV] = 1.16 x 10- 2 E^[MeV}B(Ml)[^}. (8) 

The B{E 2) values have been obtained from empirical ef¬ 
fective charges, e p = 1.5e, e„ = 0.5e, whereas the B{M1) 
values have been obtained with a four-parameter set of 
empirical g-factors, i.e. g £ = 5.586, g £ = —3.826 and 
g l p = 1, g l n = 0 (Honma et al. 2004). 

In addition, we have also included proton resonant cap¬ 
tures on the thermally excited target states. Since the 
first-excited state in 64 Ge is quite high (E x = 902 keV), 
thermal excitation can be neglected for typical X-ray 
burst temperatures. For the 6 ®As(p, 7 ) 66 Se rate, we in¬ 
cluded proton capture on the first four thermally excited 
states of 65 As (i.e., on the 0.187, 0.501, 0.863 and 0.947 
MeV states listed in Table 1). Capture on thermally ex¬ 
cited states contributes at most 38% to the total capture 
rate at 2 GK (Lam et al. 2016). The properties of 65 As 
and 66 Se for the ground-state captures are summarized 
in Table 1 and Table 2, respectively. In addition, the 
properties of 66 Se for the first-excited-state capture (the 
major thermally excited-state contribution) are summa¬ 
rized in Table 3. 

Peak temperatures in recent hydrodynamic XRB mod¬ 
els have approached 1.5-2 GK (Woosley et al. 2004; Jose 
et al. 2010). At such temperatures, resonant rates for 
the 64 Ge(p, 7 ) and 65 As(p, 7 ) reactions are expected to 
be dominated by levels with E r < 2.5 MeV (i.e., Gamow 
energy (Rolfs & Rodney 1988)). This means that exci¬ 
tation energy regions of up to E x < 2.5 MeV for 65 As, 
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TABLE 2 

Properties of 66 Se for the ground-state capture utilized in the present 65 As(p,7) 66 Se resonant rate calculation. 


■Tf M 

E x [MeV] 

E T P l b l £* heo I 3 ! 

E r [MeV] [c] 

c 2 s 7/2 
(l — 3) 

C 2 S 3/ 2 
{l = 1) 

c 2 s 5/2 

(1 = 3) 

c 2 s 1/2 
(1 = 1) 

r 7 [eV] 

r p [eV] 

6 l >7 [eV] 

of 

0.000 

0.000 



0.746 






2f 

0.929(7) 

0.982 

— 

0.037 

0.069 

0.119 

0.045 

2.765 X 10~ 4 

— 

— 

of 


1.130 

— 


0.013 



2.020 X 10 -8 

— 

— 

2f 


1.552 

— 

0.002 

0.011 

0.071 

0.051 

1.402 x 10~ 4 

— 

— 

°3 


1.575 

— 


0.002 



2.056 X 10 -6 

— 

— 

4 


1.952 

0.232 

0.001 

0.001 

0.180 

0.000 

4.040 X 10~ 4 

3.079 X 10~ 2 ° 

1.924 x 10~ 2 ° 

0+ 


1.989 

0.269 


0.002 



1.441 x 10 -8 

4.797 X 10~ 18 

5.996 X 10~ 19 

2 t 


2.053 

0.333[e] 

0.002 

0.005 

0.027 

0.025 

2.218 X 10~ 4 

3.517 X 10 -14 

2.198 X 10 -14 

4 f 

2.064(3) [d] 

2.110 

0.344 

0.008 


0.004 


8.382 x 10~ 4 

2.301 X 10 -16 

2.589 X 10~ 16 

3 f 


2.102 

0.382 

0.001 

0.000 

0.017 


3.835 X 10 -5 

1.057 X 10~ 14 

9.249 x 10~ 16 

2f 


2.277 

0.557[e] 

0.000 

0.004 

0.015 

0.001 

1.751 X 10~ 4 

2.205 X 10~ 9 

1.378 x 10~ 9 

3f 


2.419 

0.699 

0.000 

0.004 

0.063 


2.361 X 10~ 4 

1.781 X 10- 7 

1.557 x 10 -7 

V 


2.474 

0.754[e] 


0.016 

0.000 

0.052 

3.579 X 10 -3 

1.152 x 10~ 6 

4.306 x 10 -6 

4+ 


2.498 

0.778 

0.001 


0.115 


2.653 X 10~ 4 

1.888 X 10~ 7 

2.122 x 10 -7 

2^~ 

Z 6 


2.556 

0.836[e] 

0.001 

0.021 

0.033 

0.032 

9.915 X 10 -4 

6.723 X 10 -5 

3.935 X 10 -5 

Q + 

0 3 


2.668 

0.948 

0.000 

0.065 

0.145 


1.550 X 10 -4 

7.969 X 10~ 4 

1.136 X 10' 4 

4 l 


2.740 

1.020 

0.002 


0.040 


8.404 x 10 -4 

8.287 X 10” 6 

9.232 x 10 -6 



2.781 

1.061 [e] 


0.060 

0.351 

0.002 

5.338 X 10 -3 

4.483 X 10~ 3 

9.137 X 10 -4 

2^ 


2.865 

1.145 

0.000 

0.000 

0.000 

0.028 

1.175 X 10 -3 

5.315 X 10 -3 

6.014 x 10 -4 

if 


2.867 

1.147 


0.065 

0.001 

0.222 

2.170 X 10 -3 

5.887 X 10 -2 

7.847 X 10 -4 

3 f 


2.882 

1.162 

0.002 

0.003 

0.002 


1.209 X 10 -3 

9.439 X 10 -4 

4.639 X 10 -4 

of 


2.888 

1.168 


0.629 



1.716 X 10 -4 

1.985 X lO^ 1 

2.143 X 10 -5 

4 


2.907 

1.187 

0.000 


0.005 


3.139 X 10 -4 

1.001 X 10 -5 

1.091 X 10 -5 

4 


2.949 

1.229 

0.000 

0.003 

0.005 

0.016 

9.006 X 10 -4 

1.085 X 10 -2 

5.197 X 10 -4 

3f 


2.969 

1.249 

0.000 

0.006 

0.000 


7.804 x 10~ 4 

4.817 X 10 -3 

5.877 X 10~ 4 

4 


2.998 

1.278 

0.000 


0.079 


1.392 x 10 -3 

4.711 X 10~ 4 

3.960 X 10 -4 

4 


3.114 

1.394 

0.000 


0.157 


1.374 x 10 -3 

3.262 x 10 -3 

1.088 X 10 -3 

2 f 


3.217 

1.497 

0.000 

0.003 

0.021 

0.036 

1.534 x 10 -3 

3.221 x 10 _1 

9.539 X 10" 4 

of 


3.231 

1.511 


0.002 



2.561 X 10 -3 

2.244 x 10~ 2 

2.873 X 10~ 4 

3f 


3.237 

1.517 

0.000 

0.000 

0.006 


1.586 X 10 -3 

5.258 X 10 -3 

1.066 X 10 -3 

If 


3.274 

1.554 


0.005 

0.072 

0.010 

3.113 X 10 -3 

2.111 x 10 _1 

1.151 X 10 -3 

4 f 


3.299 

1.579 

0.000 


0.008 


9.407 X 10 -4 

1.075 X 10 -3 

5.644 x 10" 4 

3f 


3.328 

1.608 

0.001 

0.004 

0.022 


1.139 X 10 -3 

8.762 x 10~ 2 

9.842 x 10' 4 

2+ 


3.330 

1.610 

0.008 

0.083 

0.001 

0.012 

2.089 X 10 -3 

2.222 x 10+° 

1.304 x 10 -3 

2 fi 


3.357 

1.637 

0.001 

0.000 

0.027 

0.114 

1.348 X 10 -3 

2.761 X 10+° 

8.423 X 10~ 4 

5f 


3.394 

1.674 

0.000 




1.355 X 10 -3 

1.460 X 10~ 4 

1.812 x 10' 4 

4+ 


3.424 

1.704 

0.000 


0.018 


1.417 X 10 -3 

5.639 X 10~ 3 

1.274 x 10 -3 

o + 

0 8 


3.487 

1.767 

0.000 

0.000 

0.056 


4.881 X 10 -3 

2.792 x 10 -2 

3.635 X 10 -3 

o+ 

z 12 


3.499 

1.779 

0.006 

0.060 

0.001 

0.041 

4.335 X 10 -4 

7.440 X 10+° 

2.709 X 10~ 4 


[a]calculated by the present shell model; 

[bjmeasured by Obertelli et al. (2011) and Ruotsalainen et al. (2013); 

[cjcalculated by E r = E x — S p with S p = 1.720 MeV (AME2012); 

[djcalculated E r and T p based on the experimental value of E x = 2.064 MeV for this state; 
[ejresonances dominantly contributing to the rate within temperature region of 0.2—2 GK. 


and up to E x < 4.2 MeV for 66 Se should be consid¬ 
ered in the resonant rate calculations for 64 Ge(p, 7 ) and 
65 As(p, 7 ), respectively. In the present shell-model calcu¬ 
lations, the maximum excitation energies considered for 
65 As and 66 Se are 1.07 MeV and 3.50 MeV, respectively 
(see Tables 1 and 2). For 64 Ge(p, 7 ) only resonances up 
to £' r =1.035 MeV contribute significantly to the reaction 
rate up to 2 GK; for 65 As(p,j) only five resonances (at 
£V=0.333, 0.557, 0.754, 0.836 and 1.061 MeV) dominate 
the total resonant rate within the temperature region of 
0.2-2 GK. Those 21 resonances above £’ r =1.061 MeV 
make only negligible contributions to the total reaction 
rate up to 2 GK. Therefore, the contributions from the 
levels presented in Tables 1 and 2 should be adequate to 
account for these two resonant rates at XRB tempera¬ 
tures. 


2.2. Direct-capture rates 

The nonresonant direct-capture (DC) rate for pro¬ 
ton capture can be estimated by the following expres¬ 
sion (Angulo et al. 1999; Schatz et al. 2005), 


NaWv) 1 dc = 7.83 x 10 9 


xexp 


-4.249 


zr_ y /a 

UTi) 

(t) 


Sbc(Eo) 

1/3H 


[cm 3 s 1 mol 1 ], 

(9) 


with Zt being the atomic number of either 64 Ge or 65 As. 
The effective astrophysical fi-factor at the Gamow energy 
Eq, i.e., Sbc(E o), can be expressed by (Fowler et al. 
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TABLE 3 

Properties of 66 Se for the first-excited-state capture utilized in the present 65 As 5 / 2 i (p,7) 66 Se resonant rate. 


Ji M 

E x [MeV] 

E e x xp [b] El heo [a] 

E r [MeV] [c] 

c 2 s 7/2 

6 = 3 ) 

C 2 S%/2 
(1 = 1) 

c 2 s 5/2 

6 = 3 ) 

c 2 s 1/2 

6 = 1 ) 

r T [eV] 


r p [eV] 


[eV] 




1.952 

0.045 

0.002 

0.035 

0.000 

0.013 

4.040 

X 

10~ 4 

4.401 

X 

10-56 

1.834 

X 

10-56 

ot 


1.989 

0.082 



0.244 


1.441 

X 

10 -8 

2.398 

X 

IO -40 

1.998 

X 

10” 41 

2 1 


2.053 

0.146 

0.000 

0.019 

0.088 

0.098 

2.218 

X 

10” 4 

4.499 

X 

IQ -26 

1.875 

X 

10-26 

< 

2.064(3) [d] 

2.110 

0.157 

0.001 

0.016 

0.002 


8.382 

X 

10 -4 

1.261 

X 

10-25 

9.458 

X 

10-26 

3 ^ 


2.102 

0.195 

0.001 

0.011 

0.000 

0.001 

3.835 

X 

10 -5 

4.039 

X 

10-22 

2.356 

X 

10-22 

2 1 


2.277 

0.370 

0.001 

0.038 

0.020 

0.184 

1.751 

X 

10- 4 

4.549 

X 

10 -12 

1.895 

X 

IO - 12 

3 1 


2.419[e] 

0.512 

0.001 

0.083 

0.004 

0.226 

2.361 

X 

10 -4 

1.862 

X 

10 -8 

1.086 

X 

10 -8 



2.474 

0.567 

0.000 

0.000 

0.002 


3.579 

X 

10 -3 

5.149 

X 

IO" 11 

1.287 

X 

io- 11 

4 


2.498 

0.591 

0.000 

0.043 

0.208 


2.653 

X 

10 -4 

6.133 

X 

10 -8 

4.599 

X 

10 -8 

2+ 

z 6 


2.556 

0.649 

0.001 

0.026 

0.191 

0.008 

9.915 

X 

io- 4 

3.292 

X 

10 -7 

1.371 

X 

10 -7 

Q + 

°3 


2.668[e] 

0.761 

0.004 

0.107 

0.000 

0.002 

1.550 

X 

10 -4 

2.983 

X 

10 -5 

1.459 

X 

io -5 

4 3 


2.740 

0.833 

0.002 

0.011 

0.001 


8.404 

X 

10 -4 

1.353 

X 

io -5 

9.987 

X 

io -6 

4 


2.781 

0.874 

0.001 

0.032 

0.002 


5.338 

X 

10 -3 

9.809 

X 

io -5 

2.408 

X 

io -5 

2 ?- 


2.865[e] 

0.958 

0.000 

0.002 

0.324 

0.054 

1.175 

X 

10 -3 

6.788 

X 

IO" 4 

1.793 

X 

io- 4 

i i 


2.867 

0.960 

0.000 

0.000 

0.004 


2.170 

X 

10 -3 

2.502 

X 

io- 7 

6.254 

X 

10~ 8 

3 i 


2.882[e] 

0.975 

0.001 

0.008 

0.001 

0.035 

1.209 

X 

10 -3 

7.087 

X 

10 -4 

2.607 

X 

IO -4 

0 1 


2.888 

0.981 



0.089 


1.716 

X 

io- 4 

8.636 

X 

10 -6 

6.852 

X 

10 -7 

4 


2.907 

1.000 

0.000 

0.000 

0.518 


3.139 

X 

io- 4 

6.881 

X 

10 -5 

4.233 

X 

10 -5 

2£ 


2.949 

1.042 

0.000 

0.006 

0.012 

0.007 

9.006 

X 

io- 4 

6.740 

X 

HW 4 

1.606 

X 

10' 4 



2.969[e] 

1.062 

0.002 

0.003 

0.001 

0.025 

7.804 

X 

io- 4 

1.754 

X 

10 -3 

3.151 

X 

io- 4 

4 


2.998[e] 

1.091 

0.003 

0.025 

0.000 


1.392 

X 

10 -3 

2.795 

X 

10 -3 

6.970 

X 

io- 4 

4 


3.114[e] 

1.207 

0.001 

0.061 

0.103 


1.374 

X 

10 -3 

3.138 

X 

10 -2 

9.872 

X 

io- 4 

2+ 


3.217 

1.310 

0.000 

0.019 

0.016 

0.001 

1.534 

X 

10 -3 

3.314 

X 

10~ 2 

6.107 

X 

io- 4 

°6 


3.231 

1.324 



0.025 


2.561 

X 

10 -3 

2.545 

X 

io- 4 

1.929 

X 

10 -5 

4 


3.237 

1.330 

0.001 

0.007 

0.000 

0.053 

1.586 

X 

10 -3 

1.026 

X 

10 _1 

9.112 

X 

io- 4 

4 


3.274 

1.367 

0.002 

0.002 

0.002 


3.113 

X 

10 -3 

4.840 

X 

10 -3 

4.737 

X 

io- 4 

4 


3.299 

1.392 

0.000 

0.003 

0.026 


9.407 

X 

io- 4 

1.039 

X 

10~ 2 

6.470 

X 

io- 4 

3 i 


3.328 

1.421 

0.001 

0.073 

0.003 

0.019 

1.139 

X 

10 -3 

4.284 

X 

10 _1 

6.629 

X 

10 -4 

2 fo 


3.330 

1.423 

0.000 

0.015 

0.036 

0.001 

2.089 

X 

10 -3 

7.650 

X 

10 -2 

8.473 

X 

10 -4 

A 


3.357 

1.450 

0.000 

0.020 

0.001 

0.006 

1.348 

X 

10 -3 

1.496 

X 

10 _1 

5.568 

X 

10 -4 

5 i 


3.394 

1.487 

0.002 


0.000 


1.355 

X 

10 -3 

3.592 

X 

10+ 3 

1.242 

X 

10 -3 

4 


3.424 

1.517 

0.000 

0.013 

0.145 


1.417 

X 

10 -3 

1.591 

X 

10 _1 

1.054 

X 

10 -3 

3 i 


3.487 

1.580 

0.000 

0.051 

0.000 

0.001 

4.881 

X 

10 -3 

9.800 

X 

10 _1 

2.833 

X 

10 -3 

2 + 

"12 


3.499 

1.592 

0.001 

0.006 

0.041 

0.008 

4.335 

X 

10“ 4 

2.659 

X 

10 _1 

1.803 

X 

10 -4 


[a]calculated by the present shell model; 

[bjmeasured by Obertelli et al. (2011) and Ruotsalainen et al. (2013); 


[c] calculated by E r = E x — S p — 0.187 in units of MeV, with S p = 1.720 MeV (AME2012); 

[d] calculated E r and T p based on the experimental value of E x — 2.064 MeV for this state; 
[ejresonances dominantly contributing to the rate within temperature region of 0.2—2 GK. 


1964; Rolfs & Rodney 1988), 

Sbc(Eo) = ^(°)(l + I^)’ ( 10 ) 

where ^(O) is the S'-factor at zero energy, and the di¬ 
mensionless parameter r is given numerically by r = 
4.2487(Z|,/z/Tg) _1 ' /3 for the proton capture. 

In this work, we have calculated the DC S’-factors with 
the RADCAP code (Bertulani et al. 2003). The Woods- 
Saxon nuclear potential (central + spin orbit) and a 
Coulomb potential of uniform-charge distribution were 
utilized in the calculation. The nuclear central poten¬ 
tial Vo was determined by matching the bound-state en¬ 
ergies. The spectroscopic factors were taken from the 
shell model calculation and are listed in Table 1 and Ta¬ 
ble 2. The optical-potential parameters (Huang et al. 
2010) are Rq = R s . 0 . = Rc = 1-25 x (1 + At ) 1 / 3 fm, 
a o = a s . 0 . = 0.65 fm, with a spin-orbit potential depth 


of T4.o. = —10 MeV. Here, i?o, Rs. o., and Rc are the 
radii of central potential, the spin-orbit potential and 
the Coulomb potential, respectively; do and a s . 0 . are the 
corresponding diffuseness parameters in the central and 
spin-orbit potentials, respectively. 

For the 65 As(p, 7 ) 66 Se reaction, 5(0) values for DC cap¬ 
tures into the ground state and the first-excited state 
(E x =929 keV) in 66 Se are calculated to be 8.3 and 3.5 
MeV-b, respectively. The total DC rate for this reaction 
is only about 0.1% that of the resonant one at 0.05 GK. 
For the 64 Ge(p, 7 ) 65 As reaction, we find a DC 5(0) value 
for this reaction of about 35 MeV-b. The DC contribu¬ 
tion is only about 0.3% even at the lowest temperature 
of 0.06 GK. Even when considering estimated upper lim¬ 
its to the DC contribution (He et al. 2014), the resonant 
contributions still dominate the total rates above 0.06 
GK and 0.05 GK for 64 Ge(p, 7 ) 65 As and 65 As(p, 7 ) 66 Se 
reactions, respectively. The probabilities of populating 
the first-excited states in 64 Ge (£4=902 keV) and 65 As 
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(E x =187 keV) relative to the ground states at tempera¬ 
tures below 0.1 GK are extremely small, and hence con¬ 
tributions of the direct-capture from these excited states 
can be neglected. 

3. RESULTS 


The resulting total thermonuclear 64 Ge(p, 7 ) 65 As and 
65 As(p, 7 ) 66 Se rates are listed in Table 4 as functions of 
temperature. The present ( Present , hereafter) rates can 
be parameterized by the standard format of Rauscher & 
Thielemann (2000). For 64 Ge(p, 7 ) 65 As, we find: 


J 


13 819 12 211 

N a (<jv) = exp(-78.204---+ ' . q + 81.566T q 1/3 - 13.138T 9 + 1.717T q /3 - 16.149 InTg) 

t 9 T g 1/3 

+ exp(— 93.260 - 10,059 + 15-169 + 61.887T q 1/3 + 21.717T g - 8.625T q /3 - 22.943 InTg) 
T g t 9 1/3 

3 788 1Q 347 

+ exp(— 75.104 - —— + — tt— + 52.700T q 1/3 - 32.227Tg + 14.766T q /3 + 1.270 InTg), 
T g T q ' 3 


with a fitting error of less than 0.3% at 0.1-2 GK; for 
65 As(p, 7 ) 66 Se, we find: 


Na(ctv) =exp(-111.177- —— + 5 i J + 106.669T q 1/3 - 64.623T 9 + 13.521T q 5/3 + 31.256 InTg) 

Tg Tg /3 

+ exp(— 124.702 - 12,436 + 52,765 + 89.593T g 1/3 - 12.219Tg + 0.456T 9 5/3 - 2.886 InTg) 

Tg r 9 1/3 

r 2Q2 424 

+ exp(— 116.814— -+ q + 48.281T q 1/3 + 83.320Tg - 188.849T q /3 + 21.362 InTg), 

Tg Tt' 


with a fitting error of less than 0.4% at 0.1-2 GK. We em¬ 
phasize that the above fits are only valid with the stated 
error over the temperature range of 0.1-2 GK. Above 2 
GK, one may, for example, match our rates to statistical 
model calculations (see e.g., NACRE by Angulo et al. 
(1999)). 

Figure 1 shows the comparison of the Present 
64 Ge(p, 7 ) 65 As rate with others compiled in JINA REA- 
CLIB: rpsm , rath , thra, laur , and ths8. Note that only 
the rpsm rate uses an S p ( 65 As) value that is within ler of 
the recently determined experimental value. The Present 
rate differs significantly from others in the temperature 
region of interest in XRBs. The disagreement, in particu¬ 
lar with the rpsm rate, demonstrates that the statistical- 
model is not applicable for this reaction owing to the low 
density of excited states in 65 As. 

Similarly, the comparison of the Present 65 As(p, 7 ) 66 Se 
rate with other rates available in the JINA REACLIB: 
rpsm , rath , thra , laur, and ths8 , is presented in Fig. 2. 
Only the laur and rpsm rates use S’ p ( 66 Se) values that 
are within lcr of the currently accepted value. Although 
the Present rate differs significantly from the others, es¬ 
pecially at lower and higher temperature regions, it is 
consistent with all others within the remaining large un¬ 
certainties. At T >1 GK, the laur rate is the lowest rate 
simply because only three excited states were considered 


by Van Wormer et al. (1994). It should be noted that 



T [GK] 


Fig. 1.— Reaction rates of the ® 4 Ge(p, 7 ) 65 As reaction (in units 
of cm 3 mol -1 s —J ). The Present rate (red line) together with the 
upper and lower limits deduced from uncertainties are shown by 
a (green) colored band. Other available rates from JINA REA¬ 
CLIB (Cyburt et al. 2010) are shown for comparison. See details 
in the text and Table 4. (A color version of this figure is available 
in the online journal.) 
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TABLE 4 

Thermonuclear 64 Ge(p,7) 65 As and 65 As(p,7) 66 Se rates in units of cm 3 s _1 mol _1 . 


t 9 

N a (<?v) 

64 Ge(p, 7 ) 65 As 

lower 

upper 

N A {av) 

® 5 As(p, 7 )®®Se 

lower 

upper 

0.2 

4.79 x 10- 17 

3.99 x 10 _2 ° 

1.23 x 10“ 15 

1.86 x 10“ 16 

9.97 x lO’ 18 

1.45 x 10“ 15 

0.3 

1.33 x 10“ 13 

8.44 x lO' 15 

7.66 x 10“ 13 

2.21 x 10“ 12 

1.79 x lO’ 12 

1.92 x 10 -11 

0.4 

3.08 x lO” 11 

9.92 x 10" 13 

6.56 x 10 -11 

1.96 x 10~ 9 

8.15 x 10~ 10 

8.26 x 10' 9 

0.5 

1.89 x 10" 9 

1.98 x 10" 10 

4.30 x 10 -9 

1.51 x 10“ 7 

3.03 x 10" 8 

4.65 x 10~ 7 

0.6 

5.47 x lO" 8 

1.42 x 10 -8 

8.93 x 10 -8 

3.08 x 10“® 

3.31 x 10 -7 

7.55 x 10-® 

0.7 

6.86 x 10 -7 

2.61 x 10“ 7 

9.38 x 10 -7 

2.93 x 10“ 5 

1.73 x 10 _ ® 

5.91 x 10 -5 

0.8 

4.62 x 10 _ ® 

2.15 x 10“® 

6.52 x 10-® 

1.69 x 10 -4 

6.47 x 10-® 

2.94 x 10 -4 

0.9 

2.02 x 10“ 5 

8.05 x 10"® 

3.52 x 10“ 5 

6.86 x 10 -4 

2.18 x 10“ 5 

1.08 x 10~ 3 

1.0 

6.53 x lO" 5 

2.17 x 10“ 5 

1.39 x 10 -4 

2.15 x 10 -3 

7.04 x 10 -5 

3.19 x 10~ 3 

1.1 

1.69 x 10 -4 

4.69 x 10“ 5 

4.30 x 10“ 4 

5.55 x 10 -3 

2.12 x 10“ 4 

7.96 x 10~ 3 

1.2 

3.69 x 10" 4 

8.83 x 10" 5 

1.09 x 10 -3 

1.23 x 10 -2 

4.67 x 10 -4 

1.74 x 10 -2 

1.3 

7.10 x 10 -4 

1.50 x 10“ 4 

2.38 x 10“ 3 

2.42 x 10“ 2 

9.89 x 10 -4 

3.41 x 10 -2 

1.4 

1.24 x 10 -3 

2.34 x 10“ 4 

4.62 x 10“ 3 

4.33 x 10" 2 

1.90 x 10“ 3 

6.11 x 10" 2 

1.5 

1.98 x 10 -3 

3.42 x 10“ 4 

8.16 x 10“ 3 

7.18 x 10" 2 

3.33 x 10" 3 

1.02 x 10" 1 

1.6 

2.98 x 10 -3 

4.73 x 10“ 4 

1.33 x 10 -2 

1.12 x 10 _1 

5.69 x 10 -3 

1.59 x 10" 1 

1.7 

4.25 x 10“ 3 

6.28 x 10“ 4 

2.04 x 10“ 2 

1.65 x 10” 1 

9.38 x 10" 3 

2.34 x 10” 1 

1.8 

5.78 x 10" 3 

8.02 x 10“ 4 

2.97 x 10 -2 

2.34 x 10" 1 

1.49 x 10 -2 

3.31 x 10 _1 

1.9 

7.58 x 10" 3 

9.93 x lO" 4 

4.13 x 10 -2 

3.18 x lO” 1 

2.22 x 10 -2 

4.50 x 10 _1 

2.0 

9.63 x 10" 3 

1.20 x 10 -3 

5.52 x 10 -2 

4.19 x lO’ 1 

3.21 x 10 -2 

5.93 x lO’ 1 


the shell model calculation provides the first reliable es¬ 
timate of the uncertainty of the 65 As(p,y) reaction rate, 
especially as the Hauser-Feshbach rates may suffer from 
unknown systematic errors due to the limited applicabil¬ 
ity of the statistical model near the proton drip line. 

Uncertainties for the Present 64 Ge(p,y) 65 As and 
65 As(p,y) 66 Se rates were estimated by considering the 
uncertainties in the S p values (±85 keV for 65 As and 
±310 keV for 66 Se) and estimated uncertainties in the 
calculated level energies (±168 keV for both 65 As and 



T [GK] 


Fig. 2.— As Fig. 1, but for the ® 5 As(p, 7 )®®Se reaction. Note 
the rpsm rate is quite close to the rath , ths8 and thra rates at 
different temperature regions (Cyburt et al. 2010), and hence the 
corresponding line is not clearly visible. See details in the text and 
Table 4. (A color version of this figure is available in the online 
journal.) 


66 Se (Honma et al. 2002) 2 ). These were added in 
quadrature to give uncertainties of ±188 keV and ±353 
keV for the resonance energies E r of 64 Ge(p,y) 65 As and 
65 As(p,y) 66 Se, respectively. Note that for the two known 
levels, i.e., E x =187 keV in 65 As and E x =2064 keV in 
66 Se, an experimental excitation energy uncertainty of 
±3 keV is used instead. For 64 Ge(p,y) 65 As, all reso¬ 
nance strengths wy are proportional to T p since r p <cr 7 
(see Table 1 and Eq. 3). The uncertainties in r p (or wy) 
owing to the uncertainties in E r are calculated based on 
the energy dependence expressed in Eq. 5. In the case of 
65 As(p,y) 66 Se, only five resonances (at E r =0.333, 0.557, 
0.754, 0.836 and 1.061 MeV dominate the resonant rate 
within the temperature region of 0.2-2 GK. Here, the 
resonant strengths wy are proportional to T p for the first 
four resonances, while the strength of the last resonance 
at 1.061 MeV (dominating over ~0.9-2 GK) depends on 
both, r p and T 7 . Uncertainties in T 7 can be neglected 
because of the much larger rate uncertainties caused by 
E r and T p . The uncertainties for all resonances listed in 
Table 2 are considered in the present calculations. 

4. ASTROPHYSICAL IMPLICATION 

We examine the impact of our new 64 Ge(p,y) and 
65 As(p,y) rates and their uncertainties on the rp -process 
using one-zone XRB models. Post processing calcula¬ 
tions using temperature and density trajectories from 
the literature enable a quick assessment of the impact 
of nuclear physics changes on the strength of the 64 Ge 
waiting point using the A = 64 abundance, and on the 
burst energy generation rate. We use the post-processing 
approach for the K04 X-ray burst model (Parikh et al. 

2 The 66 Zn case is studied with the present model space and 
interaction, and an rms deviation between the experimental and 
calculated level energies is found to be about 140 keV (Lam et al. 
2016). 
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Fig. 3.— Final abundances, as mass fractions X , following 
one-zone XRB calculations using the K04 thermodynamic his¬ 
tory (Parikh et al. 2008, 2009). Results using rates determined in 
the Present work and in JIN A REACLIB: rpsm , rath , thra , laur, 
and ths8 are indicated. Panels (a) and (b) show the effects of using 
either different 64 Ge(p, 7 ) or different 65 As(p, 7 ) rates, respectively, 
while panel (c) shows the effect of using different 64 Ge(p, 7 ) and 
65 As(p, 7 ) rates together. The impact of the uncertainties in the 
Present 64 Ge(p, 7 ) and 65 As(p, 7 ) rates (see Figs. 1 and 2) is indi¬ 
cated as A (Present) in panels (a) and (b). (A color version of this 
figure is available in the online journal.) 


2008, 2009). However, postprocessing calculations do 
not take into account the changes in temperature and 
density that result from the energy generation changes. 
They can therefore not predict reliably the quantitative 
impact on produced abundances and light curves. To 
account for this effect, we also use the full one-zone X- 
ray burst model (Schatz et al. 2001), which represents a 
more extreme burst with very hydrogen rich ignition. 

4.1. Post processing results for K04 

With the representative K04 thermodynamic his¬ 
tory (Parikh et al. 2008, 2009), final abundances (as mass 
fractions X ) and the nuclear energy generation rate E gen 
during a burst have been studied by performing separate 
XRB model calculations with different rates. In the K04 


model, the peak temperature T pea k = 1.4 GK is similar 
to those reached at the base of the envelope in compara¬ 
ble hydrodynamic XRB models (e.g., 1.3 GK in Jose et 
al. (2010)). Figs. 3 and 4 compare results for X and E gen 
using rates from the present work to results using rates 
available in JIN A REACLIB: laur, rath, rpsm, thra, ths8. 
The impact of (a) using different 64 Ge(p, 7 ) rates (with 
the 65 As (p, 7 ) rate held constant at the Present value), 
(b) using different 65 As(p,'y) rates (with the 64 Ge(p, 7 ) 
rate held constant at the Present value) and (c) using 
different 64 Ge(p, 7 ) and 65 As(p, r y) rates together, is in¬ 
dicated in each of the two figures. For each change in 
reaction rate, the corresponding inverse reaction rate is 
also changed to maintain detailed balance. This inverse 
rate strongly depends on the adopted reaction S p -value 
for the respective forward rate. As we compare the im¬ 
pact of different rates that have been determined using 
very different S p values (see Sect. 2 and Figs. 1 and 2), 
the results illustrate not only the influence of the rate 
calculation, but also the influence of different forward to 
reverse rate ratios due to different S p - values. 

The results of Figs. 3 and 4 are interesting, but not 
entirely unexpected. For the case of the 64 Ge(p, 7 ) re¬ 
action, an equilibrium between the rates of the forward 
(p, 7 ) and reverse (7 ,p) processes is quickly established 
due to the small (p,j) Q-value relative to kT at XRB 
temperatures: at 1 GK, kT r; 100 keV. As a result, it 
is the (jp, 7 ) Q-value, rather than the actual 64 Ge(p, 7 ) 
rate, that is the most important nuclear physics quan¬ 
tity needed to characterize the equilibrium abundance 
of 65 As (and the subsequent flow of material to heavier 
nuclei through the 65 As(p,'y) reaction), c.f. Schatz et al. 
(1998); Iliadis (2007); Parikh et al. (2009). This is nicely 
illustrated in Fig. 3(a): the 64 Ge(p, 7 ) rates adopting pos¬ 
itive (p, 7 ) Q-values (thra, laur, rath ) give relatively lower 
final abundances around A = 64 and larger abundances 
at higher masses precisely because of the larger equi¬ 
librium abundances of 65 As during the burst, allowing 
for increased flows of abundances to higher masses via 
the 65 As(p, 7 ) reaction. On the other hand, the opposite 
is true for those rates adopting negative (p,y) Q-values 
(ths8, Present, rpsm ) because of the larger equilibrium 
abundances of 64 Ge and lower relative abundances of 
65 As during the burst. Indeed, the summed mass frac¬ 
tions of species with A > 70, X70, vary considerably for 
different choices of Q- values. For example, when the thra , 
ths8 or Present rates are adopted, X 70 = 0.58, 0.21, or 
0.33, respectively. 

A consequence of the increased flow of abundances to 
heavier nuclei is seen in Fig. 4(a), where the models 
adopting the thra, laur, and rath rates give the largest 
.Egen at late times due to energy released from the decay 
of the larger amounts of heavy nuclei produced during 
the burst. As expected, the opposite is true for E gen 
in the models using the ths8, Present and rpsm rates. 
We note that the predictions for E gen at late times vary 
rather significantly between the models using the differ¬ 
ent rates, with differences as large as a factor of ~ 2 . 

For the case of the 65 As(p,j) reaction, where a large 
positive (p, 7 ) Q- value is adopted in all rate estimates, the 
importance of the actual rate is illustrated in Fig. 3(b). 
The model adopting the largest rate at the most rele¬ 
vant temperatures (T > 0.9 GK, c.f. Jose et al. (2010)), 
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Fig. 4.— As Fig. 3, but for nuclear energy generation rates E gen during the burst in the K04 XRB model. Panels (a) and (b) show 
the effects of using either different 64 Ge(p, 7 ) or different 65 As(p, 7 ) rates, respectively, while panel (c) shows the effect of using different 
64 Ge(p, 7 ) and 65 As(p, 7 ) rates together. The impact of the uncertainties in the Present 64 Ge(p, 7 ) and 65 As(p, 7 ) rates (see Figs. 1 and 2) 
is indicated in panels (a) and (b). Panels to the right show expanded views of the panels to the left. (A color version of this figure is 
available in the online journal.) 


rath , gives the lowest abundances around A = 64 and 
the largest abundances at higher masses. Again, as ex¬ 
pected, the opposite is true for the model using the lowest 
65 As(p, 7 ) rate, laur. The variation in X^q for different 
choices of the 65 As(p,j) rate is significant, with X 70 = 
0.49 with the rath rate, and 0.29 with the laur rate. The 
behavior of E gen for these models is again in accord with 
the distributions of the final abundances, with the largest 
Eg en at late times arising from the model using the rath 
rate, and the lowest E gen arising from the model using 
the laur rate. 

Finally, the effects of using different 64 Ge(p,y) and 
65 As (p, 7 ) rates from the same theoretical model calcu¬ 
lation are shown in Figs. 3(c) and 4(c). This reveals 
the impact of competing influences from these two rates. 
For example, the model using the ths 8 64 Ge(p, 7 ) rate 
gave the lowest relative abundances at higher masses (see 
Fig. 3(a)), while the model using the ths 8 65 As(p,'y) rate 
gave among the highest relative abundances at higher 
masses (see Fig. 3(b)). When these two rates are used 
together, the combined impact on the final abundances 


(and E gen ) is, not surprisingly, moderated. 

Fig. 3(c) also shows that using the Present 64 Ge(p, 7 ) 
and 65 As(p, 7 ) rates results in the strongest 64 Ge wait¬ 
ing point, and the lowest final abundances at the high¬ 
est masses. Differences with respect to predictions using 
rates from JINA REACLIB are as large as a factor of 
~ 7 at individual values of A. X-jq calculated with the 
two Present rates differs by as much as factor of 1.8 from 
models using the other rates. 

We have also examined the impact on XRB model pre¬ 
dictions of our uncertainties in the Present 64 Ge(p, 7 ) and 
65 As(p, 7 ) rates, as shown in Figs. 1 and 2. Reverse rates 
for the lower and upper forward rates were determined 
using exactly the Q- values adopted for the corresponding 
forward rate calculations. Panels (a) and (b) of Figs. 3 
and 4 show how these rate uncertainties affect final abun¬ 
dances and -Egen in the K04 model, with mass fractions 
above A = 64 varying by up to a factor of 3 due to the 
individual uncertainties in the rates, X 70 varying by up 
to a factor of 2, and E gen varying by up to 35% at late 
times. The impact on X and E gen of the uncertainties 
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Fig. 5.— Luminosity per gram of material calculated with the 
one-zone X-ray burst model for nuclear physics input that, within 
uncertainties, maximally favours (initially high solid line) and 
maximally disfavours (initially low solid line) the rp -process flow 
through the 64 Ge waiting point. The dashed lines show the results 
when only the 65 As(p, 7 ) reaction rate is changed within uncertain¬ 
ties, and all other nuclear physics input is fixed. (A color version 
of this figure is available in the online journal.) 



Fig. 6.— Final mass fractions, summed by mass number, cal¬ 
culated with the one-zone X-ray burst model for nuclear physics 
input that, within uncertainties, maximally favours (open black 
circles, solid line) and maximally disfavours (filled red circles) the 
rp-process flow through the 64 Ge waiting point. (A color version 
of this figure is available in the online journal.) 

in the Present rates is somewhat smaller than that from 
different choices of rates but clearly not insignificant. As 
such, the mass of 66 Se should be determined experimen¬ 
tally and the uncertainty in the mass of 65 As should be 
reduced to better constrain model predictions. 

4.2. One-zone X-ray burst model results 

In addition we explored the impact of the remaining 
nuclear physics uncertainties related to the 64 Ge waiting 



Fig. 7.— Ratio of the final mass fractions shown in Fig. 6. 


point on the full one-zone X-ray burst model described 
in Schatz et al. (2001). We note that this model is dif¬ 
ferent from K04. It represents a burst that ignites in 
a very hydrogen rich environment, for example at high 
accretion rates and very low accreted metallicity, and 
was developed to explore the maximum extent of the re¬ 
process towards heavy elements. To determine the total 
remaining uncertainty in the burst model due to the nu¬ 
clear physics of the 64 Ge waiting point, we performed 
two extreme calculations. The calculations assume the 
most favourable (unfavourable) nuclear physics choice for 
the rp-process to pass through the 64 Ge waiting point, 
adopting the upper (lower) limit of the 64 Ge(p, 7 ) and 
65 As(p, 7 ) reaction rates, and the upper (lower) limits 
of S'p( 65 As) and S' p ( 66 Se). Varying S p values indepen¬ 
dently, rather than varying individual masses, is justi¬ 
fied as the uncertainties in S' p ( 65 As) and S p ( 66 Se) are 
each completely dominated by the mass uncertainty of 
65 As and 66 Se, respectively. Figs. 5 and 6 show the im¬ 
pact of the nuclear physics uncertainties on burst light 
curve and final composition. Clearly the nuclear physics 
uncertainties have a strong impact on observables. The 
abundance ratio of A = 64 to A = 68, a measure for 
the strength of the 64 Ge waiting point varies from 1.1 
(making 64 Ge the strongest waiting point) to 0.2 (mak¬ 
ing 64 Ge not a significant waiting point). This is con¬ 
sistent with the result from Tu et al. (2011) who found 
that with their new 65 As mass 64 Ge is only a weak wait¬ 
ing point. However, as we show here, when taking into 
account all nuclear physics uncertainties, a strong 64 Ge 
waiting point cannot be ruled out. Fig. 7 shows the ratio 
of the final abundances. Similar to the results obtained 
for model K04 using post-processing, a weak 64 Ge wait¬ 
ing point (low 64 Ge abundance) leads to an enhancement 
of the production of heavier elements by up to a factor of 
2. As already noted by Tu et al. (2011) the production 
of the heaviest nuclei with A > 106 is however reduced 
for a weak 64 Ge waiting point. This somewhat counter 
intuitive result is a consequence of the faster burning at 
higher temperature, which leads to a shorter burst as 
hydrogen is consumed more quickly. This effect is not 
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seen in the post-processing calculation because there the 
temperature trajectory is fixed. 

We also investigated the relative contributions of the 
various nuclear physics uncertainties. Similar to the K04 
post processing results we find that the 64 Ge(p, 7 ) re¬ 
action rate itself, for a fixed Q-value, has no influence 
on the burst model due to {p,'))-{'),p) equilibrium be¬ 
tween 64 Ge and 65 As. A calculation where only S’ p ( 65 As) 
and S’ p ( 66 Se) changes produced virtually the same result 
as the full variation, demonstrating that the mass un¬ 
certainties currently dominate. Changing each S p sepa¬ 
rately indicates that the 85 keV uncertainty of £ p ( 65 As) 
due to the 65 As mass uncertainty, and the 310 keV un¬ 
certainty of S' p ( 66 Se) mainly due to the unknown 66 Se 
mass contribute roughly equally. However, varying the 
65 As(p, 7 ) reaction rate within our new uncertainties, 
while leaving S p ( 65 As) and 5 p ( 66 Se) fixed at their nomi¬ 
nal values, still led to significant light curve changes (see 
Fig. 5) and a change of the A = 64 to A = 68 ratio from 1 
to 0.7. This shows that once the mass uncertainties are 
addressed, the 65 As(p, 7 ) reaction rate uncertainty will 
still play a role (even though fixing S’ p ( 66 Se) will reduce 
the rate uncertainty somewhat). 

5. SUMMARY AND CONCLUSION 

We have determined new thermonuclear rates for the 
64 Ge(p, 7 ) 65 As and 65 As(p, 7 ) 66 Se reactions based on 
large-scale shell model calculations and proton separa¬ 
tion energies for 65 As and 66 Se derived using measured 
masses and, for 66 Se, the AME2012 extrapolation. These 
rates differ strongly from other rates available in the lit¬ 
erature. For example, at « 1 GK, our 64 Ge(p, 7 ) rate is 
up to a factor of « 90 lower than other rates, while our 
65 As(p, 7 ) rate differs by up to a factor of « 3 from other 
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rates. 

We also determined for the first time reliable uncer¬ 
tainties for the 64 Ge(p, 7 ) 65 As and 65 As(p, 7 ) 66 Se reac¬ 
tions. We find that in two different X-ray burst models, 
the remaining uncertainties in 5' p ( 65 As), Sp^Se), and 
the 65 As(p, 7 ) reaction rate lead to large uncertainties in 
the strength of the 64 Ge waiting point in the rp-process, 
the produced amount of A = 64 material in the burst 
ashes that will ultimately decay to 64 Zn, the produced 
amount of heavier nuclei beyond A = 64, and the burst 
light curve. These effects are robust and appear in two 
different X-ray burst models. We conclude that to ad¬ 
dress these uncertainties a more precise measurement of 
the 65 As mass, a measurement of the 66 Se mass, and a 
measurement of the excitation energies of states in 66 Se 
that serve as important resonances for the 65 As(p, 7 ) 66 Se 
reaction will be important. 
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